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We present an exact quantization condition for the time independent solutions (energy eigen¬ 
states) of the one-dimensional Dirac equation with a scalar potential well characterized by only two 
‘effective’ turning points (defined by the roots of V{x) -\- me? = ±E) for a given energy E and 
satisfying mc^ + min V(x) > 0. This result generalizes the previously known non-relativistic quan¬ 
tization formula and preserves many physically desirable symmetries, besides attaining the correct 
non-relativistic limit. Numerical calculations demonstrate the utility of the formula for computing 
accurate energy eigenvalues. 


I. INTRODUCTION 

In this paper we present a quantization formula for 
the one-dimensional Dirac equation using the analyti¬ 
cal transfer matrix (ATM) method, thereby extending 
the non-relativistic analogue for the Schrodinger equa¬ 
tion found by Cao et al.^“^. We focus on the relativistic 
energy levels offered by a simple confining scalar poten¬ 
tial well with two turning points. 

The high accuracy of the eigenvalues computable from 
this formula together with other applications like ground 
state reconstruction should prove to be useful in many 
applications of the Dirac equation, especially in the con¬ 
text of solid state physics^’ 

Our result can be thought of as a completion of well 
known semi-classical quantization formulae (like Bohr- 
Sommerfeld and WKB) in the following sense. The semi- 
classical quantization formulae provide reliable estimates 
of the energy eigenvalues only in the limit of large quan¬ 
tum numbers, while the ATM-quantization formula gives 
exact eigenvalues for all quantum numbers. As a result, 
Cao et al. additionally characterize this formula as an 
‘exact quantization formula.’ 

In anticipation of the derivation given below, we would 
like to pursue this comparison a little further. The reader 
will recall that, in the WKB quantization formula, for 
instance, the boundary conditions at the turning points 
when applied to a suitably chosen approximate ansatz 
(for the actual wave function), lead to quantization of 
energy. However, the ATM method is based on transfer 
matrices which, in the appropriate limit ‘recovers’ the 
exact wave function dl and the (corresponding) energy 
eigenvalue E. Structurally, the exact quantization for¬ 
mula contains the usual WKB term with a non-trivial 
correction contributed by the so-called ‘scattered sub- 
waves’, which is of great significance (discussed in Sec¬ 
tion II). 

Particularly, for the Dirac equation the negative en¬ 
ergy (antiparticle) solutions make the energy spectrum 
unbounded from below, which for general potentials is 
difficult to account for using the ATM method. Also, the 
generalization of this method for even the Schrodinger 
equation to potentials having more than two classical 


turning points (x for which V(x) = E) is not clear at 
present. This difficulty translates to our inability to ob¬ 
tain the quantization condition for the Dirac problem for 
potential wells that either (1) give more than two classi¬ 
cal turning points, or (2) satisfy mc^ -I- min V(x) < 0, or 
both. Owing to these limitations, we narrow our focus 
to potential wells with only two classical turning points 
and satisfying mc^ + minU(a::) > 0. For this restricted 
class of potential wells we can apply the ATM method 
successfully. 

Following up with the derivation in Section II we point 
to some desirable symmetries of the quantization formula 
in Section III and discuss its non-relativistic limit. Sec¬ 
tion IV is devoted to numerical results and applications. 
We conclude in Section V outlining prospects of further 
study. 

II. FORMULATION 

Consider the one-particle Dirac Hamiltonian E[ = 
cap + f3{mc'^ + V), wherep = —ihd/dx is the momentum 
operator and m (c) is the rest mass of the particle (speed 
of light). We choose to represent the Dirac matrix cy.{f3) 
by the Pauli matrix <Ty{az), which has the advantage 

that the two-component wave function T' = ('0i ' 02 )^ 
can be chosen to be real^’®. The time independent Dirac 
equation prescribes the eigenvalue problem ElA! = EA!, 
where E is the energy of the particle. Further, the prop¬ 
erty i/'k' = —EAt' ^ At' = cTj.'I' shows that the positive 
and negative solutions occur in pairs. Hence, it suffices 
to consider positive energy solutions alone. 

The (reduced) Compton wavelength h/mc and rest 
mass energy mc^ provide natural length and energy scales 
in this problem. Thus, expressing every quantity in di¬ 
mensionless form X —> h/mc E ^ mc^e, V —>• mc^'d, 
the Dirac equation translates into the coupled first order 
system of equations 

•01 = (1-f i?(f)-f e)'02 (1) 

02 = (1 + I?(f) - e)^! (2) 

the overhead dot denoting differentiation w.r.t. We 
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FIG. 1. The (scalar) potential well showing effective 

turning points ^l,r for a typical energy e. 


identify solutions of the equations 1 + 'd(^) = ±£ as effec¬ 
tive turning points. For the schematic potential well de¬ 
picted in Fig. 1, two turning points < ^r) are ob¬ 

tained for a typical energy e. Note that more than two ef¬ 
fective turning points would result for |e| < —(1 -|-mini?) 
unless the restrictions noted in Section I are enforced on 
F. 

Next, we partition the interval between the two turning 
points into N segments < ^2 ■■■< < ^n+i = 

Denoting the width of the n**' segment by (5„, we 
choose an arbitrary point (a tag) within each segment 
and replace the scalar potential by a piecewise-constant 
approximation for which the potential in the n segment 
is given by dn = T^Un)- Consequently, on this segment 
the system of Eqs. (1) and (2) admits a general solution 
of the form 


_ A f ^ ^ piK„(C-Cn) I n f ^ ^ 

(3) 

where An and Bn are arbitrary coefficients, and 

Note that in the region of interest /v„, A„ are real, hence 
the wave function components ■i/'i 2 oscillatory. In 


favor of a simpler notation, we drop the superscript n 
from the wave function and infer the segment label from 
the context. 

Ensuring the continuity of 'F" at the ends of the seg¬ 
ment leads to the ( transfer) matrix equation 

V'l(Cn)^ ^ COs(KnSn) /l/,^(^n+l)\ 

i’2{^n) J \A„sin(K„5„) COs{KnSn) J yp2{Cn+l) J 

(4) 

Since we allow only two turning points at this stage, A„ 
is nonzero in each segment. Consequently, the transfer 
matrices are always well defined. Now, left multiplying 
Eq. (4) by {—fp2{^n) and dividing the resulting 

equation by ij^i{^n)'4’ii^n+i), we arrive at 

(Pn 1) ( ^ 

^ ^ \XnSm{Kn5n) C0s(k„(5„)/ \-Pn+l 

where Pj = —'02(Ci)/'0i('?i)- Simplifying this matrix 
equation gives the recurrence formula 

^ = tan ^ tan”^ (“7^) ~ 

^ tan“^ ~ tan“^ = ztt — 

z = 0,l,2,... (6) 

Rearranging Eq. (6) and summing over n from 1 to 
yields 



N N-1 , 

RnSn + y] tan“^ ( 

—! \An+l 

n=l n—1 ^ 


— tan ^ 


Pn+l \ 

A„ J 


= ZTT -|- tan ^ 


Pn+i \ 
Xn ) 



z = 0,l,2,... 

(7) 


The desired quantization condition emerges as a limit 
of Eq. (7) as S{:= maxi5„) —>■ 0. In the limiting event, 
the continuous potential variation i?(^) is recovered, with 
Pi —>■ P(Cl) and Pn+i P{^r). In Appendix A we 
prove that (for bound state wave functions) P{^l) < 0 < 
P(^_r) and |P(^L,fi)| < 00 . Also, as 5 0, A^at -)> 

X{^l,r) = 0. As a result, the ‘half phase losses’ at the 
effective turning points are given by 


lim tan ^ 

^Pn+A 

= — lim tan ^ I 


S->-0 ' 

V Aat y 

S-i-O ' 

+ 1 / 


S ■■= max(5„. (8) 


Next, we account for the phase contribution of the so- 
called scattered sub- waves. Defining the phase contri¬ 
bution for each segment 


A(j)n ■= tan 


P 


n+1 

^n+l 


— tan 


P 


n+1 


= tan ^ 


I Pn+l (An An+l) 

V Pn+l ffi AnAn-l-1 


Pn-H (An-l-1 An) ^ ^ 
Pn+l ffi AnAn +1 
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which results from expanding the inverse tangent in a 
Taylor series in powers of (A„+i — A„), we obtain the 
total phase contribution of the scattered sub-waves by 


N-l 


N-1 


lim > Afpn = — lim > 
5-s-O ^ ^ 


n—1 

N-l 


PA 


/ \ PA 

lim^0((A„+r-A„f)=- / (9) 

n=l 

Thus, in the limit 5 —>■ 0, Eq. (7) takes the form 




'iL 


PX 


P2 + A2 


d^ = (z-hl)7r, z = 0,1,2... (10) 


We can phrase this result differently by substituting 
2 ^ 

£ + 7?(0 + l’ 

and using Eqs. (1) and (2), obtaining 


( 11 ) 


P : _ _ V>i'i/>2(1 + ^(C) + g) _ 

+ A2 (1-1- -d{^) + £)^| - (1 + '&{£,) - £)'01 

^ ^ £^(C) 

A(l+l?(0+g)^ At \4>lil2 — i’24’1 J 

em ( 

2At (^«'t(-icrj^d/d0^y ^ ’ 

_ £At ^ 

4i\/£^ - At^ V'1'1 (-*'’'y'l/d^) ^a J 

■: 2icr^ = [ay,(T;^]. (13) 


Using this result and replacing the Pauli matrices with 
the Dirac matrices gives a representation independent 
form (Ref. Section Ill) of Eq. (10). Eurther, restoring 
the dimensions of the physical quantities and using n 
(instead of z) for the quantum number yields 



dx J \'lAt(cap)dA J 

= (n + l)7r, 


n = 0,l,2,... 


(14) 


where K = y/E"^ — (mc^ + V{x))^, and the effective 
turning points satisfy V(xl,r) + = ±E. 


III. DISCUSSION 

We begin by describing two symmetries of Eq. (14) 
that are physically desirable. 


A. Representation independence 

In Section II we chose a convenient representation of 
the Dirac matrices ol{( 3) — (Ty{az). However, the quan¬ 
tized energy levels are a property of the potential U(x), 


hence should be independent of the chosen representa¬ 
tion. This feature is already built into the quantiza¬ 
tion condition and can be shown in the following way. 
We know that the Dirac matrices satisfy the algebra: 

= 0^ = TL and {a, f3} = 0. Given any other represen¬ 
tation of the Dirac matrices (denoted a',0) satisfying 
the same algebra, (a generalization of) Pauli’s fundamen¬ 
tal theorem® asserts that there exists a unique (invert¬ 
ible) matrix S (up to a multiplicative complex constant) 
such that 


7 = 5-S'5, 7 = a,/3 (15) 

A direct substitution shows that this transformation pre¬ 
serves the structure of the Dirac equation (for the same 
E) with the wave function transforming as 4 a' = 54a. Ad¬ 
ditionally, as the Dirac matrices are hermitian, we must 
have 5 I = S~^ (unitarity of S). Using the transfor¬ 
mation Eq. (15) in Eq. (14) and noting that \a,0\ = 
5l[a',/3']5 gives 

4't[a,/3]4' 4''l[a',/3']4'' 

4At (cap) 4 a 4A't(ca'p)4A' 

thus confirming the representation independence of the 
obtained quantization condition. 


B. Symmetry w.r.t. the sign of E 

Although we assumed E to be positive in the derivation 
of Eq. (14), the quantization condition should not depend 
on the sign of the eigenvalue. Since A is a function of 
the apparent asymmetry stems from the factor 

4Atcr2,4A 

4At (—ACTj,d/d.£) 4 a 

in Eq. (12). Using the properties cr^ = era, = 
crj,, {<jx,cry} = 0 we obtain 

4'llcra;4' _ £(cra,^)lcra;(cra,4') 

^4'tl (-ACTyd/d^) 4' (cra,4')t {-icr^cTyd/d^) 4' 

_ (-£)((Ta,4')lcra;(o-a,4') 

(cr a, 4A)t {-icTyd/dC) {(Txd>) 

which in the light of our earlier observation 774' = £'4' 
i7((T a; 4a) = —£(<7 a; 4a) establishes the expected symme¬ 
try. As the effective turning points remain unchanged 
under £ >->• —£, Eq. (14) also holds for the antiparticle 
solutions. 

Einally, we look at the non-relativistic limit of Eq. (14) 
and affirm that it reproduces the Schrodinger quantiza¬ 
tion formula found by Cao et al.^“® in this limit. The non- 
relativistic limit con cerns energies £ ~ mc^\ and is easily 
demonstrated by the replacements: (1) mc^-|-U(a;)-|-£ « 
2m(? and (2) me? + V{x) — E k. V{x) — Eg, where 
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Eg := E — mc^ (the ‘shifted’ energy) is a small quan¬ 
tity (Cf. Sec. 4.4 of Ref. [7])'^ . Using this prescription, 
Eq. (1) and (2) (after replacing the dimensions) become 


d^i 

dx 

d-i/^2 

dx 


2mc 


V{x) - Eg 
he 


V’l 


(16) 

(17) 


which further decouple to yield the one-dimensional 
Schrddinger equation 


dx^ 


Orri 

— {V{x)-Eg)i;, 


(18) 


From Eq. (11) we arrive at the limiting forms of A and P 


where xl^{xrJ denotes the left (right) ‘classical’ turning 
point given by the smaller (larger) of the two roots of the 
equation V (x) = Eg. The above equation is the correct 
non-relativistic limit of Eq. (14), which agrees with the 
well-known quantization formula describing the bound 
states of the Schrddinger equation^^^. 

IV. NUMERICAL RESULTS 


We turn now to a discussion of computational aspects 
of Eq. (14) and outline some applications of the same. 
Before proceeding, we point to an apparent impediment, 
namely, the integrand singularities at the turning points 
xl,r arising as a result of 


E dK 
- E2 dx 


.E dV 
^ K dx’ 


K{xl,r) = 0. 


( 22 ) 


A 


Eg - V{x) 
2mc^ 


P'. 


2mc 


-Ps 


(19) 


with 


1 d^i 
ijji dx 


( 20 ) 


where the subscript s is intended for notational homo¬ 
geneity of later equations. Proceeding further, we sub¬ 
stitute these results into Eq. (10) making the following 
observations 


(Eg - V{x)) 


me 

h 12m 


me 


n? 


{Eg - V{x)) := - Kg 


P dA 


Pg 


2^] 

+ A2 de '’ p 2 + ^(^Eg- V{x)) \ dx) 


Pg 


dU/da 




+ 1^1 \ V2mc2 {Eg - V{x)) 


4z\h^iEg-V{x)) 


Pg + Kg \ dx V mc2 

h Pg dK 


me Pg + Kg dx 


thus arriving at (d^ = me/h dx) 
dKg / Pg 


dKg 

dx 




dx « ( 2 -I-l)7r, z = 0,1,2... 

( 21 ) 


Without further qualification, such points might cause 
the integral to diverge. Luckily, this does not happen, as 
we rigorously show in Appendix B. Furthermore, conver¬ 
gence requirements do not restrict the class of admissi¬ 
ble potentials to which our analysis applies. With that 
caveat, we look at the problem of determining energy 
eigenvalues of the Dirac equation for a given scalar po¬ 
tential U(x). 

The collection of one-dimensional potentials for which 
the Dirac equation is closed-form solvable is rather small. 
Even for the modest (scalar) ‘simple’ harmonic oscillator, 
the wave-functions cannot be obtained in closed analytic 
form. To our knowledge the one-dimensional Woods- 
Saxon potential®, the linear confining potential®’® and the 
scalar exponential potential^® have enjoyed closed form 
solutions so far. Thus, an accurate quantization formula 
becomes particularly useful. 

Consider the linear confining potential V = g\x\, g > 0 
whose wave functions are written in terms of the Hermite 
function H^{x)^^. Adopting the conventions h = e = 1, 
V = E"^ I2g given in Ref. [9], we wish to incorporate 
the known wave functions into Eq. (14) and reproduce 
the energy eigenvalues, which are otherwise obtained 
as solutions of the transcendental equation H^{a) = 
2i/E[^_^{a), where a = mf^/g. 

Note that the authors of®’® solve this eigenvalue prob¬ 
lem in the so called Jackiw-Rebbi representation of the 
Dirac equation corresponding to a = cr^ and /3 = ax in 
which case Eq. (14) takes the form 



K 


EV 


V’N V’i \ 
' 01'02 - ' 0201 / 


dx = {n + l)7r. 


(23) 


where the overhead dot now denotes differentiation w.r.t. 
X. Substituting the wave functions in the l.h.s of Eq. (23), 
we obtain 


r.i /A-7 . . -i( « \ f Hl{x)-2vHl_^{x) \ 

Ia{v) =-Kv — ayJ2v — — 2i'tan , = + / -twtt -rw , - 

\\/2v - ) Ja \H^{x) - H^-i{x)Hi,+i{x) J ^2v - x'^ 


dx 
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Ia{v) 

TT 



FIG. 2. Plot of for a = 1 (continuous line ) a = 2 

(broken line). 

the derivation of which makes use of the recurrence for¬ 
mula = 2xHi,{x) — 2i'H^-i{xy^. The quanti¬ 

zation condition can be rewritten as Ia{v) = {n + l)7r. 
Note that for E < mc^ + min V there can be no turn¬ 
ing points. Therefore, bound state solutions can only be 
expected for z/ > a^/2. 

In Figure 2 we plot the graphs of /a(z/)/7r for a = 1, 2 
as a function of ly. These curves intersect the horizon¬ 
tal integer- lines (shown in the same hgure) precisely at 
the locations of the energy eigenvalues, thus confirming 
remarkably the validity of Eq. (14). For comparison, we 
collect the graphical intersection points (computed using 
bisection search) in Table I where the exact eigenvalues 
obtained by Cavalcanti® are also listed. 

Now, concerning potentials for which the wave func¬ 
tions are not obtainable in terms of known functions, a 
numerical routine can be outlined for the Dirac prob¬ 
lem (following the prescription given in^^ for the non- 
relativistic ATM-quantization formula) yielding eigen¬ 
value estimates, whose accuracy is only appreciable for 
larger eigenvalues-a feature that we do not understand 
completely at present. Therefore, a discussion of this 
scheme is not given here. 

Next, we briefly mention another interesting applica¬ 
tion, to motivate further work. In fact, a detailed ac¬ 
count of the same will be addressed in a follow- up paper. 
Particularly, we consider the problem of reconstructing 
the ground state wave function Iko, for a given potential 
V(x). The basic scheme is to assume the ground state 
energy E^, to be a parameter and invert Eq. (14) to ob¬ 
tain 'Po(Eo). Although this is easier said than done, for 
simple power-law potentials V(x) = g\x\'^{g > 0, n G N), 
using the ansatz 


ing powers of Eq on either side of the resulting equation. 
We would then resort to a variational principle to mini¬ 
mize the energy functional E[4'], which is the expectation 
of the Dirac Hamiltonian w.r.t. the obtained paramet¬ 
ric wave function, to find Eq- Similarly, recovery of the 
potential with knowledge of the ground state wave func¬ 
tion is conceivable. However, this inversion problem is 
significantly complicated by the presence of the effective 
turning points, which depend on the potential implicitly. 

V. CONCLUSION 

In this paper we analyzed the bound states of the one¬ 
dimensional Dirac equation for a scalar potential well 
V(x) (satisfying two constraints laid in Section I), obtain¬ 
ing an exact energy quantization formula that extends 
the previously known analog for the Schrddinger equa¬ 
tion. In fact, we could show that our formula reproduces 
this result in the appropriate non-relativistic limit. Fur¬ 
ther, we discussed the physically desirable symmetries of 
the quantization formula and applied the same to com¬ 
pute the energy eigenvalues for the potential V = g|a;|. 

Before concluding, we invite interested readers to strive 
for a completion of our relativistic quantization formula 
that would apply to an arbitrary scalar potential. Such 
a generalization must address two cases for which our 
analysis fails. First, the case of non-confining poten¬ 
tials: consider, for example, the one-dimensional Dirac 
equation with the scalar exponential potential V(x) = 
Ae~^^, A > 0^® which, despite being repulsive every¬ 
where, supports bound states for A < 0. This is a conse¬ 
quence of the nature of coupling in the Dirac equation. 
Note that the exponential potential yields only one effec¬ 
tive turning point, thus failing to satisfy ( at least) one of 
the constraints laid in Section I. As a result, our analysis 
would not apply to this case. Second, one must account 
for potentials that yield more than two effective turning 
points (a double well potential for instance). Unfortu¬ 
nately, the generalization of the analytic transfer matrix 
method to such potentials even for the non-relativistic 
problem stands open. 
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P = exp 




in Eq. (10), one can obtain the coefficients aj by compar- 


Appendix A 

Let T be a bound state wave function, with its com¬ 
ponents ipi ^2 satisfying Eqs. (1) and (2). Based on the 
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TABLE I. First five values of ly = /2g for select values of a = m/^fg. {h = c = 1). 


V 

a = 1 


Q = 2 


Exact 

ATM 

Exact 

ATM 

VQ 

1.39627444057259’* 

1.39627444809303 

3.338595 40177509’* 

3.33859536647797 


3.05676024015993 

3.05676024192944 

5.452160 76495126 

5.45216075601056 

v-i 

4.30627664769999 

4.30627665789798 

7.006087 30469830 

7.00608729608357 

V-i 

5.61521082352847 

5.61521084997803 

8.568945 86286508 

8.56894588172436 

1^4 

6.80477121323347 

6.80477123537566 

9.978608 33615064 

9.97860836439766 


® Underlined digits are not given in®, and can be obtained by solving H^{a) = 


properties of an admissible bound state wave function we 
deduce an useful property of the auxiliary function 

p(e) = -^ (Ai) 

Wi 

which is well defined (and bounded) at any finite ^ ex¬ 
cepting the nodes of t/^i. We show that 

P{^l) < 0 < P(efl), (A2) 

where denotes the left (right) effective turning 

point (defined by 'd{^L,R) + 1 = e^)- Using this result 
we obtain the ‘ half phase losses’ at these turning points 
to be 7r/2 (Section II) . We recall that the chosen rep¬ 
resentation of the Dirac matrices allows us to work with 
a real Hence, the inequality in proposition (A2) is 
valid. In proving the above proposition, we require the 
following lemma. 

Lemma 1. The components o/d' cannot vanish simulta¬ 
neously at any finite Equivalently, $ ^ & R|'I'(^) = 0. 

Proof. Suppose there exists a such that d'(^o) = 0. 
Consider a finite interval containing in which the wave 
function can be represented as {A, B arbitrary coeffi¬ 
cients) 



where tp^'^ are the linearly independent solutions of 
the second order ODE 


the Wronskian of the two functions. But this would con¬ 
tradict the linear independence of the solutions 
at ^o, hence such a does not exist. □ 

The regions where d(^) -|- 1 < (>)e are designated as 
allowed (forbidden) regions. Consider the properties 

PI: V' 1,2 —t 0 as 1^1 —oo 

P2: ipi ^2 ^ 0 for any ^ (no nodes) in the forbidden region 

which hold for any bound state wave function. We now 
prove inequality (A2). 

Proof. Since Eqs. (1) and (2) remain form invariant 
under the transformation ^ i—>■ —^(=> t Cr) E^nd 

tpi 1 —^ —4>i while P I—>■ —P, it suffices to prove any 
one of the two inequalities in (A2). We focus on the 
right effective turning point The truth of proposition 
(A2) rejects the possibility sgn['0i(^/{)] = sgn[i/' 2 (Cfl)] 
where sgn[ ] is the signum function. To prove this we 
let iPi^ 2 {^r) > 0. From P2, it follows that tpi > 0 
for all ^ > ^R {a, forbidden region). Since tp 2 {^R) > 0, 
Eq. (1) implies that ■i/'i(Ci?) > 0. Since ipi is increasing at 
it must attain at least one maximum before vanish¬ 
ing asymptotically as oo (property PI), remaining 
positive-definite all along. Clearly, at the site of this 
maximum, '0i = 0 i /’2 = 0, which contradicts property 

P2 for ip 2 - The other possibility ■i/'i,2(Cfl) < 0 (a ‘reflec¬ 
tion’ of the previous case) is readily contradicted from 
form-invariance of Eqs. (1) and (2) under the transfor¬ 
mation T- -T-. Thus, sgn[V'i(Cfl)] = -sgn[V' 2 (■?/?)] =» 
P(^fl) >0. □ 


'ipi 


( m \ 

ye+m + ^j 


ipip 


f e + m + A 
\e-m-A 


tpi =0 


(A4) 


Finally, we show that 


\P{Am)\ < oo 


(A6) 


that results from decoupling Eqs. (1) and (2).’ 
Since 'k vanishes at ^o, the matrix equation 


(tp[^\^o) tp^\^o)\ 



(A5) 


must hold. Thus, in order to prevent the wave function 
from vanishing identically (i.e. A = H = 0) on the inter¬ 
val, we must have W{ip[^\^o), tp^\^o)) = 0 where W is 


Proof. Since ipi and '02 cannot vanish simultaneously (us¬ 
ing the above Lemma), we need only show that an effec¬ 
tive turning point cannot be a node of ipi. Consider 
as before. Assume tpi{^R) = 0. From PI, '0i —>• 0 
as ^ oo. Thus, ipi must attain at least one minimum 
(maximum) if' 02 (Cfl) < 0(> 0) (•.• sguji/ji] = sgn[V' 2 ]), re¬ 
maining negative (positive) definite all along. Similar to 
our earlier argument, at the site of this extremum 02 = 0, 
which contradicts property P2 . Thus, 0i {£,r) 0 0. □ 
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Appendix B 


combining these results, we obtain 


In this appendix we show that the integral (Eq. (10)) 

PA 

'-I (Bi) 


irK r_H_ 

"4 ^£2 _ (1+^)2 |P| 


exists for all potentials 'd(^) that satisfy the constraints 
laid in Section I namely, (i) mini? + 1 > 0 and (ii) 1 + 
i?(^) = £ holds only at the two effective turning points 
^L.R- the analysis that follows, it is important to bear 
in mind that the components 'tpi ,2 of a bound state wave 
function 'h have only finite number of zeros (or nodes) 
for < ^ < Furthermore, between two consecutive 
nodes of ipi there exists exactly one node of 1^2 and vice- 
versa, which is a consequence of the fact that between 
two consecutive nodes of ipi (say) there exists at least 
one extreme point of ' 0 i at which "02 oc 'i/^i = 0 (see 
Eq. (2)). Since this argument holds good for both '0i and 
^p 2 , there could at the most be one node of 'ip 2 between 
two consecutive nodes of ' 0 i. 

Now, we show that / does not fail to exist due to 
the singularity at ^r. To do so, we focus on an inter¬ 
val , ^r) where is greater than the largest among 
the nodes of ■^i, 'il ’2 and the (only) point of minimum of 
(denoted by ^c)- Defining 

riR p\ 


we have 


|/a|< 



PA 


r^R 

A 

L 

P 2 -hA 2 

de< 

L 

P 




(B3) 


Using Eq. (11) we find 


|A| 


£ _ 1^1 

1 + i ^ + e ^£ 2 - (^1 + ^)2 


Also, as 1 -I- mini? > 0 (constraint (i)). 


since i?(^) > 0 for ^ > ^c- Integrating the r.h.s. of (B4) 
by parts we obtain 


sin-1 

ij I <__f z 

' ^'-2\P{U)\ |P(^a)| 



_d 

df 



d^ 


Note that |7^(Ca)| is non-zero and finite, owing to the 
choice of the point and so is |T’(^i^)| (inequalities (A2) 
and (A 6 )). The r.h.s above can be bounded from above 
by its absolute value leading to 


141 < 


TT 

2 


1 

Wi^\ 


1 

WM 




d? 


(B5) 

Note that, sgn(P(^)) = sgn(P(^fi)) for all < ^ < Cfl, 
as P doesn’t cross the ^ axis (no nodes of i/' 2 ). Thus, 
|P| = P (see inequality (A2)). Furthermore, the deriva¬ 
tive of 1/P has no zeros in this interval, since 


d 

de 



0 ^ p2 


1 -|-1? — e 
l-hi?-h 


which refutes the fact that P is real in the chosen repre¬ 
sentation of the Dirac matrices. Therefore, 



_d 

de 




1 1 


Pi^R) Pi^a) 


Substituting this result in inequality (B5) we obtain 


141 < 


Pi^n) 


(B 6 ) 


e 

1 -I- 1 ? -I- e 


< 1 


which is finite (due to inequalities (A2) and (A 6 )). A 
similar argument works for the left turning point. 
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